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Magnetorheological (MR) dampers enable wide-range adjustment of damping and are con¬ 
sidered as promising actuators for suspension performance optimization. However, MR 
dampers also exhibit complex nonlinearities, including hysteresis and bi-viscosity. 
Moreover, in addition to handling road irregularities, the suspension system needs to 
adjust its performance according to variations in vehicle load. The above nonlinearities 
and requirements both bring challenges to the design of semi-active suspension 
controllers. In this paper, the nonlinearities of MR dampers are handled with a piecewise 
modeling method, and the vibration attenuation is realized by designing a static 
piecewise-affine (PWA) controller. Furthermore, to deal with the time-varying load 
and to fully utilize the capability of MR dampers, a parameter-dependent piecewise- 
quadratic Lyapunov function (PDPQLF) is employed, which results in a load-dependent 
PWA (LDPWA) Hoo controller. Simulation and experimental results under random, bumpy 
and sinusoidal excitations show that the proposed controller can achieve enhancements on 
both comfort and handling performance. Above all, the proposed method enables effective 
high-frequency vibration suppression for underload conditions, while for overload 
conditions, the low-frequency performance is more significantly improved. 

© 2019 Elsevier Ltd. All rights reserved. 


1. Introduction 

As key components of vehicles, suspension systems have been studied for many years, and numerous methods related to 
structures [1], actuators [2] and control [3] have been proposed to improve suspension performance. When faced with road 
surface irregularities, vehicle suspensions are required to (1) enhance the comfort performance, (2) ensure the handling per¬ 
formance, and (3) limit the suspension deflection. These requirements, however, are conflicting to each other and a tradeoff 
has to be involved [4]. More complicated, the requirements of suspensions will change with the variation of vehicle loads [5]. 
For example, a vehicle which is configured for heavy load should have a relatively stiff suspension. However, when the load is 
removed, the suspension will be excessively stiff or over-damped, which results in deteriorative comfort performance. Con¬ 
versely, when the suspension is configured for an unloaded condition, it will be too soft to support additional weight. There¬ 
fore, it is imperative to design a suspension system which can achieve the multi-objective performance for both road 
excitations and load variations. Since traditional passive suspensions cannot be adjusted online, semi-active and active sus¬ 
pensions have attracted more and more attention of industries and academia to deal with the above problems. 
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Active suspensions, which can inject energy to suspension systems by adjusting the active force directly, have powerful 
adjustment capabilities and can effectively handle the complex road and load conditions [6-8]. Besides many literatures on 
the attenuation of road disturbances, the researches of active suspension on load adaptation can be traced back to the 1990s. 
For example, by combining an air spring and an air responsive shock absorber (stiff for a high air pressure which is caused by 
a heavy load and soft for a low air pressure which corresponds to a light load), an active suspension which can respond to 
both road and load conditions is designed and patented by the Monroe corp. [9]. Further, as a modification of this concept, a 
booster which enables the pressure communication between the air spring and the shock absorber is proposed by the Ten- 
neco Inc. [10]. In addition to these suspension structure innovations, some researchers also achieve the load adaption prop¬ 
erty by designing control algorithms of the active force. For example, a load-dependent controller is creatively given in [11 ] 
where the gain matrix can vary according to the vehicle load to achieve a better multi-objective performance. This 
parameter-dependent control strategy is further expanded in [12] to design the vehicle active suspension controllers which 
consider both the vehicle inertial parameter variations and the time delays of electro-hydraulic actuators. Recently, an adap¬ 
tive controller for active seat suspensions has been proposed in [13] and the controller is experimentally verified to show 
good performance on vibration control and high system robustness to the operator’s mass variations. Above all, the load- 
adaptive (or load-dependent) active suspensions have been heavily studied in the past two decades. However, the high 
structure complexity and the huge power consumption will inevitably increase the cost and limit the application of active 
suspensions. 

Conversely, semi-active suspensions, which can only dissipate energy from the vibration system by adjusting the damp¬ 
ing force of shock absorbers, have a compact size and require very little energy supply [14-16]. Although semi-active sus¬ 
pensions are not as adjustable as active suspensions, they are more convenient for practical applications. For the above 
reasons, the development of load-adaptive semi-active suspensions has attracted the interest of industry. For example, in 
[17], an integrated damper which includes a damper and a valve block is designed. The valve block has a spring pilot which 
can indicate the pressure on the spring. This patented damper enables damping response to load conditions without relying 
on conventional air (or hydraulic) springs. However, the assembly and manufacture of this damper are not a short-term job. 
Therefore, design of a load-adaptive algorithm for semi-active suspensions has become a research trend. Currently, magne- 
torhelogical (MR) dampers which are treated as promising semi-active actuators due to their wide adjustment range and fast 
response speed are drawing more and more attention [18,19]. Moreover, this actuator can be easily driven by a current of 
0 A-3 A, which is convenient to be integrated into traditional vehicle electrical systems. However, MR dampers also show 
up complex nonlinearities including dissipative constraint, bi-viscous and hysteresis [20-22], and the existing load- 
adaptive controllers for active suspensions cannot be directly transplanted into the MR semi-active suspensions. Thus, these 
nonlinearities pose challenges to the controller design. As a pioneer, a load-dependent LPV/H 2 controller is designed for the 
MR semi-active suspensions in [23], where the time-varying load is treated as a scheduling parameter and the actuator non- 
linearities are handled by an inverse MR model. The design approaches based on inverse MR models always involve a two- 
step technique. In the first step, a damping force controller is designed based on a linear suspension model, and in the second 
step, an inverse model which reflects the dynamics of MR dampers is used to transform the damping force to the driving 
current. However, as illustrated in [24,25], this method may lead to unpredictable results and ensure neither the system sta¬ 
bility nor the suspension performance, since the capabilities of MR dampers or the constraints on damping force have not 
been concerned in the first step. For this reason, in [26], the authors of this paper have proposed a piecewise MR suspension 
model to handle the actuator nonlinearities, where the nonlinear force constraint is approximated with constant constraints 
of subsystems. Then, a piecewise affine H^ (PWA-H^) controller is designed by using a global quadratic Lyapunov function 
(GQLF). The method is concise and is convenient to be applied online. However, the employment of GQLF may produce con¬ 
servative results, since a unique Lyapunov function is used for all the piecewise subsystems. Another problem is that the 
proposed PWA-H^ controller cannot adapt to the load variation. 

In this paper, a parameter-dependent piecewise quadratic Lyapunov function (PDPQLF) is employed to design a 
load-dependent PWA-H^ (LDPWA-HoJ controller. Considering the capability of MR dampers, a piecewise constrained MR 
suspension model is firstly introduced to design the damping force controller. An inverse MR damper model is then utilized 
to generate the driving current. To handle the road irregularities, the force tracking error of the inverse model and the load 
variation, a load-adaptive robust controller is further designed by using a PDPQLF. Finally, the proposed strategy of this paper 
is verified by simulations and experiments under different external excitations. The main contributions in this paper are 
summarized as follows. 

• Through the piecewise modeling of MR semi-active suspensions, the difficulty of concerning the realizability of damping 
force in the classic inverse model based control design can be solved. Additionally, compared with other models with 
nonlinear constraints [27,28], the proposed model is more convenient for the off-line design and the on-line application 
of suspension controllers. 

• With the employment of PDPQLF which empowers the controller a load-adaptive ability, the damping force in each sub¬ 
system can act up to its own constraint condition to realize the optimal suspension performance. Compared with tradi¬ 
tional Hqo methods which consider the unconstrained cases 29] or possess a high-order dynamic feedback controller [30], 
the proposed method can guarantee the nonlinear damping force constraint with a simple static feedback controller. 
Additionally, the proposed method can utilize the capabilities of MR dampers more comprehensively than the method 
in [26]. 
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• The effectiveness and applicability of the proposed controller in this paper are verified through a hardware implementa¬ 
tion in accord with actual vehicle electronic systems and bench tests under sinusoidal excitations. 

The remainder of this paper is organized as follows. Section 2 will introduce the piecewise model of semi-active suspen¬ 
sions and the control problem. Furthermore, in Section 3, a detailed design process of the LDPWA-Hoo controller will be given. 
Simulation and experimental results will be discussed in Sections 4 and 5, respectively. 

2. Problem statement 

2.1. Piecewise model of semi-active suspensions 

The vertical vehicle dynamics, especially the two dominant modes, sprung mass bouncing and wheel hopping, can be rep¬ 
resented by the quarter of vehicle (QoV) model shown in Fig. 1(a), where m s is the sprung mass which represents the vehicle 
chassis; m u is the unsprung mass which represents the integrated mass of wheel hub, tire, MR damper, spring and linkages; 
k 5 is the spring stiffness; c s is the damping of an additional passive damper; F is the damping force output by the MR damper; 
k t stands for compressibility of the pneumatic tire; z s and z u is the vertical displacements of sprung and unsprung masses, 
respectively; and z r is the excitation when a vehicle runs on uneven road. Based on the above definitions, dynamics of MR 
semi-active suspensions can be described as follows: 

r m s z s = -k s {z s - Zu) - c s (z s - z u ) - F 
1 m u z u = k s (z s - z u ) + c s (z s - z u ) - k t (z u - z r ) + F. 

As illustrated in [9], load conditions can be divided into two types: static load and dynamics load. Static load is associated 
with the passengers of the automobile, the cargo in the automobile, and so forth. In contrast, dynamic load involves varia¬ 
tions related to road conditions and vehicle maneuvers. Similar to the researches in [11,12,23], this paper only considers the 
static load variation, which makes m s in (1) a time-varying parameter. Further, as many methods have been developed to 
estimate the vehicle load [31,32], it is supposed that m s can be measured online and resides in an interval, thus: 

( 2 ) 

To describe the dynamics on F, the velocity/force (V/F) character of MR dampers is given in Fig. 1(b) where z def = z s - z u , 
and the bi-viscosity and hysteresis are displayed. These complicate nonlinearities will inevitably increase the difficulty of 
suspension controller design. Therefore, an inverse MR model based control schema shown in Fig. 2(a) is always adopted 
in industrial applications and literatures. Accordingly, a damping force controller will give desired damping force, and the 
inverse model will change the force into the drive current of MR dampers. However, if the desired force is beyond the capa¬ 
bility of MR dampers, the design performance will be unrealizable. For this reason, a piecewise constrained suspension 
model is established in [26] to describe the force boundary. Specifically, due to the bi-viscosity, MR dampers show different 
properties at different deflection velocity levels, thus the whole force available scope can be divided into three regions, 
namely the pre-yield region and the post-yield regions shown in Fig. 1(b). In Fig. 2(b), the three yield regions are further 
divided into m ps ti, m pst 2 and m pre strip-type partitions with velocity bounds p u and p vi , and the force boundaries can be 
approximated by constant values q u and q vi , where ieC is the index of subsystems, £ := {1,..., N} and 
N = m pst i + m pst 2 + tn pre denotes the whole partition number. More details about the division and approximation can be 
referred in 26]. Accordingly, the available region of F can be declared by the following piecewise constraints: 



Fig. 1 . (a) MR semi-active suspension model; (b) Nonlinear actuator character. 
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Deflection velocity: i def (m/s) 


(a) 


(b) 


Fig. 2. (a) Inverse MR model based suspension control system; (b) Piecewise approximation of damping force boundary. 


<?Li < F < Qui, if Zdef e lPu,Puih ' e £• 


( 3 ) 


Remark 1. The velocity intervals of partitions corresponding to post-yield and pre-yield regions are suggested to choose 
discriminatively. Specifically, for the pre-yield region with a stiff damping property, a smaller interval (or a bigger m pre ) 
should be chosen to achieve a higher approximation precision. While for the post-yield regions with soft damping properties, 
bigger intervals are preferred to decrease the complexity of the piecewise model. 

Based on the above approximation, a piecewise constrained model of semi-active suspensions is derived. We define 
x = [xi x 2 x 3 x 4 ] T as the system state where X\ = z s - z u ,x 2 = z u - z r ,x 3 = z St and x 4 = z u . Then, define z = [z\ z 2 z 3 ] T 
as the performance output, where Z\ = p{i s ,z 2 = p 2 (z s -z u ),z 3 = p 3 (z u -z r ), and p x ,p 2 and p 3 are weighting coefficients. 
Furthermore, by defining u = F as the control input and w = z r as the disturbance input, a state space model can be organized 
as follows: 

x(t) = Aj(m s )x(t) + B,(m s )u(f) + W,w(f), 

z(f) = Q(m s )x(t) + Di(m s )u(t), q Li < u(t) ^ q m , x(t) e 7%, i e f 


where 7:= {x|x 3 -x 4 e [p L i,Pui]} denotes the definition region of the i-th subsystem. According to (1), the matrices in (4) 
can be obtained as follows: 



■ 0 

0 

1 

-1 ■ 



' 0 ' 


0 

0 

0 

1 



0 

A,- (m s ) = 

_ ks_ 

0 

_ _£s_ 

Cs_ 

, B,-(m s ) 

== 

1 


m s 

m s 

m s 



m s 


ks 

_ kt_ 

Cs_ 

_ _£s_ 



J_ 


. m u 

m u 

m u 

m u _ 



. _ 

D ,(m s ) = 

P\ 

m s 

0 0 

1 T 

] , W, 

= [0-10 

0] 1 



f>2 0 

- 0 p3 


P i£s 
m s 


0 

0 


P\Cs 

m s 


0 

0 


Remark 2. As the piecewise approximation influences only the control constraint, the system matrices 
(Aj(m s ),Bj(m s ), W 2 -,Ci(m s ) and D f (m s )) of different subsystems remain unchanged. 


Remark 3. In model (4), the disturbance input is completely unknown. Thus, besides the road excitation, w can also repre¬ 
sent the approximation error of the force boundary and the force tracking error. Additionally, as the force tracking is gener¬ 
ally realized through an open-loop schema (other implementations which rely on force feedback can also be referred in 
[33,34]), the inverse model accuracy will seriously influence the realization of control performance. Since the development 
of MR damper model or its inverse model is still a research frontier [35], designing a damping force controller that is robust 
to w is of practical significance. 
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Remark 4. We define x(0) = 0 when u = 0 as the system origin which is definitely the point where z def = 0 and F = 0 in 
Figs. 1(b) and 2(b). Therefore, the system partitions can be divided into two classes, thus ( = ( 0 u £i, where Co denotes the 
index set which contains the origin-included partition and Ci is the index set of origin-excluded partitions. According to 
Fig. 2(b), as the damping properties in the first and third quadrants of the V/F coordinate system are generally symmetrical, 
we can specify that q Li = -q ui for i e Co with an acceptable approximation error. 

2.2. Load-dependent piecewise control problem 

Based the above approximation, the capability of MR dampers has been quantitative described. Due to the piecewise con¬ 
stant form of the constraints, the force controller can be designed without using the MPC (model predictive control) methods 
which are usually employed to handle nonlinear constraints and show poor real-time performance [27,28]. Further, consid¬ 
ering the robust requirement mentioned above, the mature control techniques are more appropriate and practical. How¬ 
ever, the constraints of origin-excluded partitions in system (4), which are asymmetric about zero (thus q u ^ - q vi when 
i e Ci ), are difficult to be handled within the framework of H^ control design. Therefore, a PWA controller is suggested in 
[26] as follows: 

u = K*x + dj, ieC (5) 

where the affine item is predesigned as that d, = (q u + q ui )/2 to transform the asymmetric constraints to symmetric ones as 
follows: 

|K;X| d maX) i (6) 

where u max ,; = {q m - q u )/2. The feedback gain K, plays a leading role in disturbance attenuation, constraint assurance and 
system stability. 

Based on the above innovations on suspension model and controller style, the following control problem has been con¬ 
cerned in [26]. 

Problem 1. For a constant m s , define a global quadratic Lyapunov function (GQLF) as that V g (x) = x T Px. Then find a set of 
control gains K f and a matrix P, such that the closed-loop system consisting of (4) and (5) is asymptotically stable, the H 00 
norm from w to z is less than a given scalar y, and the constraints in (6) are satisfied when w T (t)w(t)dt < w max . 

The resulted GQLF-based PWA-H^ controller can guarantee the designed performance within the ability of MR dampers 
and the uniform P will lead to a concise design process. However, by depicting the design schema of Problem 1 in Fig. 3(a), 
we can see that: (1) the GQLF ensures not only the state stabilities of possible switches (adjacent partitions) but also those of 
the impossible switches (nonadjacent partitions), and (2) the GQLF is employed to handle different load conditions and dif¬ 
ferent constraints of every subsystem. Both of these aspects will lead to a conservative K, which will limit the performance of 
MR dampers. As a contrast, Fig. 3(b) shows an ideal design schema where every partition has its own Lyapunov matrix and 
the time-varying property of the Lyapunov matrices will make control gains change with the load variation. To realize this 
schema, this paper considers a PDPQLF as follows 

V p (x,m s ) =x T P,(m s )x. (7) 

Therefore, as a modification of Eq. (5), a following LDPWA controller will be designed: 

u = Ki(m s )x + di. (8) 

Continuous-time PWA systems declare that Lyapunov functions should be continuously differentiable [36], which is a 
strong requirement and is adverse to the conservatism reduction. Therefore, this paper will consider discrete-time systems. 
By choosing a suitable sampling interval, the closed-loop system consisting of (4) and (8) can be discretized as follows: 


P 




(a) 


(b) 


Fig. 3. State transition and design schemas: (a) constant GQLF-based control; (b) time-varying PQLF-based control. 
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x(k + 1) = 0j(m s )x(/<) + W iW(k) + B ,-(m s )di 
z (k) = O i{m s )x{k) + D i(m s )di 


where 0,(m s ) = A,(m s ) + B I (m s )Ki(m s ), 0,-(m s ) = Cj(m s ) + D I (m s )Ki(m s ). The system matrices A(m s ),B(ms),C(m s ),D(m s ), and 
W are the discrete-time forms of A I (ms),B ! (m s ),Ci(m s ),D ! (m s ), and W„ respectively. Further, it is obvious that the load- 
dependent system matrices of the i-th subsystem in (9) satisfy that 


A i(m s ) B i(m s ) 
Q(m s ) D i(m s ) 


2 


ZSn,( m s) 

"i=1 


Ci,„, 


D(,m 


( 10 ) 


where <5i(m s ) -f S 2 (m s ) = 1;A,^,B l; i,Q.i and D*i are constant matrices with m s = mi;Ai )2 ,Bj )2 ,Ci )2 and D i2 are those with 
m s = m 2 . Moreover, with the measurable m s ,S-i(m s ) and S 2 (m s ) can be obtained as that: 


<5i(m s ) 


mi(m 2 -m s ) d (m) = nhjnhjzUhl 

m s (m 2 -m 1 )’ ^ m s (m 2 - m,)' 


( 11 ) 


Therefore, the control problem concerned in this paper can be drawn as follows: 


Problem 2. For time-varying m s , find a set of load-dependent Lyapunov matrices Pj (m s ) for (7) and control gains K,-(m s ) for 
(8), such that system (9) is asymptotically stable with the given H^ index y and the constraints in (6) are satisfied when w is 
bounded by w max . 

To facilitate the problem solving, some definitions and conditions are pre-specified. Specifically, when the states of sys¬ 
tem (9) transit from 7 to 7 Zj at time k, we specify that the system is governed by the dynamics of the subsystem in 7Z t at that 
time. Let Q represent the set of index pairs which denote the possible switches from one partition to itself or another par¬ 
tition, that is: 

Q= {(i,j)|x(k) g TZi,x(k+ 1) g 1Zj,i g (}. (12) 

Further, the state regions TZi of origin-excluded partitions can be exactly described by following degenerate ellipsoids: 

£> = {x(k)|||F,x(k) -f/j|| < 1}, i g Ci (13) 

where F, = 2e/(p M - p Li ),f, = -(p w + p u )/(p ui - p u ), and e = [0 0 1 - 1], 


3. Controller design 

In this section, we focus on solving the problems derived in last section with linear matrix inequalities (LMIs). 

3.2. Design of the GQLF-based PWA-H^ Controller 

For Problem 1, the constant control gains in (5) can be obtained through the following theorem. 

Theorem 1. An admissible set of control gains K* in (5) exists such that the requirements in Problem 1 are satisfied, if there exist 
matrices Q 1 = Q, > 0, Y; > 0, and scalars < 0, such that the following LMIs hold: 


Qd + DfYi 



i £ Co ^ 


Hj + >^Bjdj(B;dj) ★ 

QQ, + D/Y, -I 

fiBjdjfJ + QFjf 0 



<0, ieCu 


(14a) 


(14b) 


7 2 u max,i 

W max 





(14c) 


where E* = QAj + A,-Q, + Y^B^ + B ; Y, + y 2 W ; W[, A ( , B,. Q, D, are constant system matrices, I is an identity matrix and ★ denotes 
the items induced by matrix symmetry. Moreover, the constant gains and the global Lyapunov matrix can be obtained as that: 

K,- = Y/Q 1 , P = <T\ ie£. (15) 


Proof. Derivation of this theorem has been given in [26]. □ 
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3.2. Design of the PDPQLF-based LDPWA-H^ Controller 

For Problem 2, the time-varying control gains in (8) can be obtained through the following theorem. 

Theorem 2. For all (i, j) e Q, n { e {1,2} and nj e {1, 2}, an admissible set of gains K* (m s ) in (8) exists such that the requirements 
in Problem 2 are satisfied , if there exist matrices: Qln, = Q-i.n, > 0, Qjn, = Qj.n, > 0, Gj n ., R i n . and scalars: rj i > 0, ^ > 0, such that 
the following LMIs hold: 


" —y 2 I 

★ 

★ 

★ 





Wi 

-Qjnj 

★ 

★ 








< i £ C(b 


(16a) 

0 

0 

-I 

★ 





0 

'I'n,, 







' -y 2 1 
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W ; 
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-I 


★ ★ 
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m 

(16b) 
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di) T 
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1 

cs 

1 

* 
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n ini 





u 2 








max,i 

y 2 Wmax 

Yi.n, * 

s: 0, 

i 6 Co, 



(17a) 

R„, -I 







u 2 








max,z 

y 2 Wmax 

Y f , nj ★ 


★ 





R, 

i„ -I 


★ 


<0, ie Cl 


(17b) 

FjGjiH 0 

-ftl 

'fi- 

>). 





where 

^!>, — (Ai jn . GiOj- + B ! n .R !n . j , 
n Ui — (Q^jG^n,- + Di n .R !n . j , 

Yi,nj = di n, ~ Gj ,rii ~ G i n ., 

— ~dj,rij ~ ^iBi,n i di(Bj j n I .di) . 

Moreover, the load-dependent control gains and Lyapunov matrices can be obtained as that: 

K 1 (m s ) = ^ ni (m s )R l - ni G^, (18a) 

n,=l 

Pi(m s ) = ^ B| (m s )drJl, ie£. (18b) 

n i=l 

Proof. As the item D I (m s )d I in (9) affects neither the system stability nor the performance, a sufficient condition for the 
satisfaction of the requirements on system stability and disturbance attenuation in Problem 2 is that the following inequality 
holds: 


V p (x(fe +1)) - V p (x(k)) + z T (k)z(k) - y 2 w T (k)w(k) < 0 (19) 

where z (k) = z(/<) - D ! (m s )d l . To facilitate the derivation, the time instant k and the time-varying parameter m s will be omit¬ 
ted in the following. 

Substitute (9) into (19), the following inequality can be obtained: 


[0/X + WjW + B/dj] T Pj [0/X + W iW + B,-di] - x T Pj-x + x T 0>j r O ! -x - y 2 w T w < 0. 


( 20 ) 
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Inequality (20) can be rewritten in the following matrix form: 


"x" 

7 

"o," 


' ©, " 


w 

1 

( 

0 

0 

w+ 

W, 

B f d, 

Pj(*) + 


—P, ★ ★ 

0 -y 2 l ★ 
0 0 0 


7 

"x" 


w 

. 

1 


<0. 


(21) 


For i e 7, the partition information expressed in (13) can be induced into (21) by using the S-procedure. Further, by using 
Schur complement twice for (21), we can get the following inequality: 


( 22 ) 


*i - AjF^Fj 
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★ 


0 

—y 2 l 
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where 2* is a positive scalar. Pre- and post-multiplying (22) by 
TO I 0 0 01 


°i = 


0 0 0 Pj 0 
0 0 0 0 1 
1 0 0 0 0 
0 0 10 0 


and oj respectively, we obtain that: 


' -y 1 1 

★ 

★ 

★ 


★ 

PjW, 

-Pj 

★ 

★ 


★ 

0 

0 

-I 

★ 


★ 

0 

0[PJ 

4>! 

P, A,F 

|Fi 

★ 

0 

(W) T Pj 

0 


-*(/; - 

ig the Schur complement again, inequality (23 
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As that (1 -f?) = 1 +/i(l -f?) fv expression (24) can be equally written as that: 
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where (p^ = P, - X t 1 (B,dj) P ( . Using Schur complement for (25), we can obtain that: 
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and post-multiplying (26) by 
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and i/fj respectively, and then applying Q, = P, 7 Qj = P, 1 and rj i = 1 , we can obtain that: 
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Then, according to Lemma 1 given in the Appendix, we can obtain that a sufficient condition that guarantees (27) is that 
there exists matrix Gj(m s ) such that the following inequality holds: 
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where 

% = (AjG,- + B,R i ) T . n, = (CjGj + D,R,) T 
r, = Q, - G, - Gj, % = -Q.j - t/M (Bid,) 1 , 


and R, = K,G,. Therefore, in the following, we just need to prove that (16b) can satisfy (28). 

Consider a fixed pair ( ij ), firstly, for each n i} multiply inequalities (16b) labeled with rij with S Uj and sum. Then, multiply 
the inequalities labeled with n t with S ni and sum. Due to that 5\ + S 2 = 1, we can finally obtain (28) once the matrices G- f and 
R, satisfy that: 


G/(m s ) — ^^n i (tn s )Gi j „., R;(m s ) — ^ (tn s )Ri,n, 


n,=l 


and K,,P/ satisfy (18). Similarly, for the origin-included partition (i e Co), we can obtain (16a) without considering the ellip¬ 
soidal state information. 

In the following, the LMIs that guarantee the control constraints will be derived. By summing (19) from 0 to 00 , we can 
obtain that: 


V P (x(k)) + XZl|z T (fc)z(fc)|| < y 2 J2w\k)w{k). 

k =0 

o the bounded distu 

Vp(x(k)) < y 2 w max . 


k=0 k= 0 

Due to the bounded disturbance energy, we can further obtain from (29) that: 

,, 2 . 


(29) 


(30) 


Further, for the origin-excluded partitions, (30) and (13) can be induced into (6) by using the S-procedure to obtain a suffi¬ 
cient condition which guarantees the control constraints as that: there exist positive scalars a* and ft such that the following 
inequality holds: 
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"fJf,- ★ 

X 

J — G(-i (Vp (x) — y 2 W m ax) — ft 

.1. 


J,F, f>- 1. 
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By assigning a* as u^ axi /(y 2 w max ) and using the Schur complement, (31) can be reformed as: 
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Pre- and post-multiplying (32) with ft = 
can obtain that: 
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Due to that Q, - G, - GT ^ -GJq, 1 G i , a sufficient condition that guarantee (33) is that: 


< 0. (34) 


Further, by using the same method which derives (28) from (16b), we can prove that (17b) guarantees (34). In a similar way, 
inequality (17aa) can be obtained for the origin-included partition without considering the ellipsoidal state information. The 
proof is completed. □ 

Remark 5. For Theorem 1, the number of LMIs that should be optimized together is 2N, while for Theorem 2, the optimiza¬ 
tion complexity is related to both the partition number and the element number in £1 For the division method given in 
Section 2, the number of switch pair is 3 N - 2. Further, due to the load variation, the number of possible switches that cor¬ 
responds to each (ij) will be n,nj. Therefore, the number of LMIs expressed in (16) and (17) that should be solved together is 
2n I n J (3N-2). 


max,i y _ 
y 2 w max 1 

★ 

★ 

R, 

-I 

★ 

F,G, 

0 



Remark 6. Besides the application in semi-active suspension control, Theorem 2 is applicable to general piecewise linear 
systems with asymmetrical constraints and time-varying parameters. 


4. Simulation analysis 

4.1. Simulation environment and control algorithms 

The parameters of the QoV model are given in Table 1. Additionally, three different sprung mass values mi,rh and m 2 , 
which respectively correspond to minimal, nominal and maximal load conditions, are chosen. 

4.1.1. MR damper models 

Various models had been proposed to represent the complex nonlinearities of the MR damper, for example, Bingham 
model [37], Bouc-Wen model [38], modified Bouc-Wen model [39], ANFIS model [40] and so on. The Bingham model is sim¬ 
ple, but the hysteresis cannot be adequately described. Therefore, this model is only employed in situations when there is a 
significant need for a simple model. The Bouc-Wen model or its modified form can accurately represent the nonlinearities of 
MR dampers. However, the model complexity is unavoidably increased with an extended number of model parameters 
which may impose difficulties in the identification. Additionally, Bouc-Wen models can only display the normal character 
of MR dampers (from z def and I to F), which will inevitably lead to the additional modeling work on the inverse character 
(from F and z def to I). The ANFIS model can accurately describe the normal and inverse characteristics, but requires a large 
amount of test data. Considering the above aspects, a hyperbolic tangent function based damper model from [41] is 
employed in this paper with the following expression: 

F = CoZdef + fcoZdef + y mr /tanh(c,z de f + klZdef) (35) 

where / is the electric current input, y mr is a proportional rate of damping force per ampere over the deflection velocity, c 0 is a 
linear viscous damping coefficient, k 0 is a linear stiffness coefficient, C\ and k i are hysteresis coefficients related to the deflec¬ 
tion displacement and velocity, respectively. Since every parameter has a specific physical meaning, expression (35) is a phe¬ 
nomenological model. Accordingly, the mechanical constitution of this model is depicted in Fig. 4(a). From (35), we can see 
that there are only five parameters to be identified, and the model structure is relatively simple such that an inverse damper 
model can be easily obtained as follows: 

/ = (F - CoZdef - koZdef)/ (y mr tanh(ciZdef + fciz de f)). (36) 

Models (35) and (36) can clearly describe the normal and inverse characters of MR dampers respectively, notably the hys¬ 
teresis and bilinear nonlinearities [42]. However, the model accuracy is limited, especially for the complex condition with a 
wide range of vibration frequency. Nevertheless, considering the trade-off between accuracy and complexity, model (35) is 


Table 1 

QoV model parameters. 


m s 



m u 

k s 

c s 

k t 

m-i 

m 

m 2 

751<g 

27kN/m 

400Ns/m 

400kN/m 

350 kg 

450 kg 

550 kg 
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(a) (b) 

Fig. 4. (a) Mechanical structure of hyperbolic tangent function based MR models; (b) Application schemas of damping force controllers. 


still widely used for the control design purpose [30], and will be used in this paper for the evaluation of damping force con¬ 
trollers. If rich damper test data can be obtained, and the computing power of hardware control unit is strong enough, the 
inverse ANFIS model can also be utilized to achieve higher force tracking accuracy. In this paper, the parameters in (35) are 
identified based on the test data of a real MR damper from the FAW-Tokico corporation. The resulted parameter values are 
that: 


c 0 = 810.78, Ci =5.78, k 0 = 620.79, k i = 20, y mr = 1200, 
and the electric current range is that I e [0 A,3 A]. 

4A.2. Control algorithms 

To get insight into the performance of the LDPWA-Hoo algorithm, this paper will simultaneously present the simulations 
of the PWA-Hoo algorithm and a classical semi-active algorithm namely SH-ADD. Application schemas of the three algo¬ 
rithms are shown in Fig. 4(b) where the plant is simulated by the interaction of expressions (1) and (35). Further, some pref¬ 
erences about the three algorithms are given below. 

(1 ) PWA-Hoo/LDPWA-Hoo! According to the identified MR damper model and the approximation method given in [26], the 
velocity bounds and the force constraints of each partition for the design of PWA and LDPWA H^ algorithms are listed in 
Table 2. In this paper, the partition number of N = 9 with m pre = m ps ti = m pst 2 = 3 is considered (Detailed analysis on the 
influence of the partition number can be referred in [26]). Further, the control gains of these two algorithms are obtained 
by solving the LMI optimization problems given in Theorems 1 and 2, respectively. In this paper, the LMI problems are han¬ 
dled by a commercial solver named MOSEK [43]. To compare the two algorithms fairly, a common set of weighting coeffi¬ 
cients are chosen as that p 1 =10,p 2 =0.1 and p 3 = 10. As PWA-H^ algorithm can only concern the constant load, the 
nominal load value will be used for the controller design. While for LDPWA-Hoo algorithm, as shown in Fig. 4(b), we assume 
that an accurate load estimation module is employed, thus m s = m s . Additionally, the parameters w max and y, which should 
be pre-specified for the design, will be treated as adjustable parameters in this paper. 

(2) SH-ADD: This is an effective comfort-oriented algorithm for semi-active suspensions [44]. Compared with the famous 
skyhook method, SH-ADD algorithm can reduce the vertical body acceleration in a wider frequency range. The algorithm can 
be expressed in a bang-bang form as follows: 

c = | Cmin, if [(z s - KZ 2 S ) St 0flZ s (Z s - KZ U ) > 0] U [(z s - KZ 2 5 ) < 0flZ s (Z s - KZ U ) > 0] 

1 c m ax, otherwise 


Table 2 

Parameters of system piecewise approximation. 


Velocity bounds (unit: m/s) 

Pu 

Pu. 

Pl3 

PL4 

Pl5 

Pl6 

PL7 

Pl8 

PL9 

-1.60 

-1.20 

—0.80 

-0.40 

-0.13 

0.13 

0.40 

0.80 

1.20 

Pui 

PU2 

PU3 

PlJ4 

Pus 

P U6 

P U7 

PU8 

PU9 

-1.20 

—0.80 

—0.40 

—0.13 

0.13 

0.40 

0.80 

1.20 

1.60 

Force bounds (unit: 1<N) 








<7li 

PL2 

PL3 

QL4 

QL5 

QL6 

QL7 

Q L8 

Q L9 

-4.88 

—4.54 

-4.19 

-3.84 

-0.50 

-1.07 

0.32 

0.66 

0.98 

<7ui 

Pu2 

QlJ3 

Q U4 

PU5 

PU6 

Q U7 

QU8 

^9 

-0.95 

—0.62 

—0.26 

1.14 

0.50 

3.96 

4.28 

4.59 

4.92 
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where c is the damping coefficient and k is a design parameter which is set as 121. As SH-ADD algorithm concerns none of 
the actuator nonlinearity, a linear relationship between c and the electric current is adopted in Fig. 4, thus c min —> I min and 
c m ax —> /max- In this paper, SH-ADD algorithm will act as a representation of bang-bang controllers. 

4.2. Conservatism analysis and parameter selection 

In this subsection, the conservatism of PWA and LDPWA algorithms will be evaluated from two aspects: (1) feasibility 
of the corresponding LMI optimization problems; (2) performance of the closed-loop systems. Since both two aspects are 
highly related to the choice of parameters w max and y, a principle for the parameter selection will also be derived. 

4.2A. Feasibility of LMI optimization problems 

The feasibility of the LMI optimization problems given by Theorems 1 and 2 will be tested under different values of w max 
and y. In this paper, the ranges of w max and y are chosen as (0.0001, 10000) and (0, 80), respectively. Then, by setting the 
intervals as that log 10 w max = 0.1 and y = 1, the feasibility results of 6400 parameter pairs defined as (w max ,y) are summa¬ 
rized in Fig. 5. Obviously, the feasible parameter range of LDPWA-H^ algorithm is bigger than that of PWA-H^ algorithm. 
Specifically, for PWA-H^ algorithm, a small y or a big w max both lead to unfeasible results. However, for LDPWA-Hoo algo¬ 
rithm, the corresponding LMI problem is feasible under almost all the tested parameters except that y is too small. Therefore, 
the above results also indicate that the complex form of Theorem 2 does not affect its feasibility. In addition, the resulted 
LDPWA-Hoo algorithm can theoretically handle more severe disturbance. 

4.2.2. Performance of closed-loop systems 

Furthermore, the influence of parameters on suspension performance will be analyzed through closed-loop simulations. 
As a smaller y indicates better disturbance attenuation performance, only the feasible parameter pairs with the minimal y for 
each w max will be evaluated here. In Fig. 5, these pairs make up the sets T\ and r 2 for PWA and LDPWA algorithms, 
respectively. Thus, each w max in T, corresponds to a control gain. To concentrate on the conservatism resulted from different 
design processes rather than outer conditions, a nominal load condition is simulated here and the vehicle is supposed to be 
running on a class C road which is suggested by the international organization for standardization (ISO) with a speed of 
50 km/h. The generation of this road profile can be referred in [45]. After running the simulation under different control 
gains, the root-mean-square (RMS) values of response signals z s ,z u -z r and z s -z u , which indicate the comfort, handling 
and suspension deflection performance, respectively, are summarized in Fig. 6. Six typical test points Al, Bl, Cl, A2, B2 
and C2 shown in Fig. 6 are selected with the damping forces displayed in Fig. 7. 

In Fig. 6, the performance resulted from PWA and LDPWA H^ algorithms will reach the same level when w max increases to 
a certain value. Specifically, the RMS values at points Cl and C2 are almost equal. This is because the w max is much higher 
than the actual disturbance energy. In other words, controllers designed to handle serve disturbances actually deal with rel¬ 
atively weak disturbances. As a result, in Fig. 7, the damping forces at Cl and C2 are both seriously compressed, which limits 
the performance of MR dampers. When w max falls to the values at Bl and B2, the comfort performance under both algorithms 
is improved. This is because that the performance of MR dampers is released, and the range of damping force at Bl and B2 is 
significantly expanded compared to that at Cl and C2. Meanwhile, the expanded damping force will inevitably make the sus¬ 
pension “softer”, which deteriorates the handling performance and extends the suspension deflection. Comparing the damp¬ 
ing forces at Bl and B2, it can be found that the force operating range of LDPWA-Hoo algorithm is significantly larger than that 



logio W max 


Fig. 5. Feasibility of PWA-H^ and LDPWA-H^ algorithms under different design parameters. 
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logio^max 

Fig. 6. Suspension performance resulted from different values of w max in set and r 2 . 





i de f (m/s) i de f (m/s) z def (m/s) 





idef (m/s) idef (m/s) i de f (m/s) 

Fig. 7. Damping force resulted from different values of parameter w max : PWA-H TO algorithm (Al, Bl, Cl); LDPWA-H^ algorithm (A2, B2, C2). 


of PWA-Hoo algorithm. Actually, the output of LDPWA-H^ algorithm occupies almost the entire available force region. There¬ 
fore, with a slight deterioration in handling performance, LDPWA-H^ algorithm achieves significantly better comfort perfor¬ 
mance than PWA-Hoo algorithm. Although the suspension deflection at B2 is bigger than that at Bl, it is still within 
acceptable limits. When w max is lowered to the values at Al and A2, the damper performance is further released. However, 
the output of PWA-Hoo algorithm in Fig. 7 exceeds the constraint (especially the output at post-yield regions), and results in a 
fairly “stiff’ suspension characteristic. Therefore, the comfort performance given by PWA-Hoo algorithm becomes very poor, 
and the handling performance is unexpectedly good. Conversely, the output of LDPWA-Hoo algorithm still satisfies the tor¬ 
tuous constraint and fills the entire reachable region. This is exactly the effect of the piecewise Lyapunov function. In sum¬ 
mary, compared to PWA-Hoo algorithm, the proposed LDPWA-Hoo algorithm can better utilize damper performance, and can 
maintain the designed performance over a larger parameter range. 

According to above analysis, the principles to choose appropriate design parameters can be derived as that: (1) the non¬ 
linear constraint of the damping force should be fully satisfied; and (2) the ability of MR dampers should be used as much as 
possible. Following these principles, the parameters of the two H^ algorithms in this paper are selected as follows: 

• PWA-H^: (parameters at point Bl) y = 32,w max = 0.05; 

• LDPWA-Hoo: (parameters at point A2) y = 5, w max = 0.001. 
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By designing the controller according to the above parameters and performing closed-loop simulations, the suspension 
responses shown in Fig. 8 can be obtained. According to the RMS and peak-to-peak (PP) values in Fig. 8(a), we can find that 
LDPWA-Hoo algorithm enhances the comfort performance more significantly than PWA-Hoo algorithm and passive suspen¬ 
sions (In this paper, passive suspensions are represented by MR suspensions driven by constant electric currents, and 0 A 
and 3 A denote the soft and hard damping properties, respectively.). The same conclusion can also be drawn from the 
time-domain results shown in Fig. 8(b). In summary, the proposed LDPWA-Hoo algorithm provides less conservatism than 
PWA-H^ algorithm, and the potential of MR dampers has been further explored. 


4.3. Analysis of load-adaptive performance 

Based on the above selected parameters, this section will analyze the load-adaptive performance of the proposed 
LDPWA-Hoo algorithm. To this end, two load conditions are considered here: 

• an underload condition with m s = m\ ; 

• an overload condition with m s = m 2 . 

For LDPWA-Hoo algorithm, we consider two situations for the above conditions: 

• {variable gain-VG): the time-varying sprung mass is accurately estimated, thus m s = m s ; 

• (constant gain-CG): the sprung mass is estimated incorrectly and the estimated value is a constant nominal value, thus 
m s = m. 


The performance of passive suspensions will also be analyzed to verify the necessity and effectiveness of load-adaptive 
properties. To evaluate this property within a concerned frequency range of 0 Hz-25 Hz, the simulations are set to drive the 
vehicle on a Class C road at a speed of 50 km/h, and the comfort and handling performance will be characterized by the PSD 
(power spectral density) results of signals z s and z u -z r , respectively. 

Fig. 9 gives the PSD results of passive suspensions. To facilitate analysis and descriptions, the entire frequency range in 
Figs. 9 is divided into a low-frequency band (LFB: 0 Hz-2 Hz), a medium-frequency band (MFB: 2 Hz-10 Hz) and a high- 
frequency band (HFB: 10 Hz-25 Hz). The MFB can also be regarded as the sensitive frequency band of the human body which 
is generally 4 Hz-8 Hz. Additionally, the resonant frequencies of sprung and unsprung masses are labeled SRF and URF, 
respectively. By comparing the resonant frequencies of soft and hard suspensions in Fig. 9, it is obviously that increase in 
load makes SRF shift to a smaller frequency point, while the URF is not influenced. In addition, according to Fig. 9(a), heavier 
load makes the comfort performance in MFB and HFB better. These phenomena are reasonable since the inertia of the sprung 
mass increases with load. Accordingly, conclusions can be drawn as that: (1) For underload conditions, it is more necessary 
to suppress vibrations in the MFB and HFB; (2) For overload condition, the performance within the LFB is more urgent. 

Figs. 10 and 11 show the PSD results of LDPWA-Hoo controlled suspensions and the passive ones, where Fig. 10 is obtained 
in the underload condition and Fig. 11 displays the overload condition. The forces output by LDPWA-Hoo controllers are given 
in Fig. 12. 


RMS and Peak values 


RMS (zdef * 2) Peak (z s ) 




(a) 


(b) 


Fig. 8. Performance of PWA and LDPWA algorithms under the chosen design parameters: (a) RMS and peak values; (b) time-domain responses. 
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LFB' ' [ ' MFB | HFB 1 



Frequency (Hz) 



(a) 


(b) 


Fig. 9. PSD responses of passive suspensions under different load conditions: (a) comfort performance; (b) handling performance. 




(a) (b) 

Fig. 10. PSD results of LDPWA-H^ controlled semi-active suspensions in underload condition: (a) comfort performance; (b) handling performance. 




(a) (b) 

Fig. 11. PSD results of LDPWA-H^ controlled semi-active suspensions in overload condition: (a) comfort performance; (b) handling performance. 


By comparing the results under the underload condition in Fig. 10, it can be found that the LDPWA-Hoo algorithm with 
variable gain has better comfort and handling performance in the MFB than the constant gain algorithm. This phenomenon 
can be explained by referring to Fig. 9(a), where the variable gain controller gives a softer damping characteristic and the 
constant gain controller exhibits steeper V/F characteristic. Since the performance within the MFB is more necessary under 
underload condition, the variable gain controller is preferred. Conversely, for the overload case in Fig. 11, it can be seen that 
the variable gain controller performs better in LFB and HFB, while the constant gain controller has a slight advantage in MFB. 
According to Fig. 12(b), this phenomenon can be explained by the fact that the damping characteristic given by the variable 
gain controller changes to stiff. This is significantly good, not only for vibration attenuation in LFB and HFB, but also for the 
improvement of cargo capacity. Based on the above analysis, the proposed LDPWA-Hoo algorithm can adapt to the load vari¬ 
ation, and exhibit desired suspension performance. 

4.4. Comparisons with PWA-H^ and SH-ADD algorithms 

In this subsection, the proposed LDPWA-Hoo controller will be further analyzed by comparing with PWA-Hoo and SH-ADD 
controllers. Since the PWA-Hoo and the SH-ADD algorithms do not have load-adaptive features, only the simulation under 
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( a ) 


(b) 


Fig. 12. Output damping force of LDPWA-H^ controllers: (a) underload condition; (b) overload condition. 


nominal load will be considered here. Both random and bumpy road profiles will be employed to test the semi-active 
suspensions. 

4.4.1. Random road simulations 

For random road simulations, vehicles equipped with semi-active suspensions or passive suspensions are designated to 
travel on Class C roads at a speed of 50 km/h. Fig. 13(a) and (b) show the PSD results for response signals z s and z u -z r , 
respectively. The time-domain responses of signal z s are given in Fig. 14(a), and (b) gives the corresponding RMS and PP val¬ 
ues. Similar results for signal z u - z r are displayed in Fig. 15. 

Comparing the PSD results of the two H x controllers in Fig. 13, it is obvious that the proposed LDPWA-Hoo algorithm 
shows better comfort and handling performance in LFB and MFB. However, in HFB, the performance of the two algorithms 



10 ° 

Frequency (Hz) 


10 1 



(a) 


(b) 


Fig. 13. PSD results of semi-active suspensions controlled by different algorithms: (a) comfort performance; (b) handling performance. 



(a) 


(b) 


Fig. 14. Comfort performance of semi-active suspensions controlled by different algorithms: (a) time-domain responses; (b) RMS and PP values. 
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(a) 


(b) 


Fig. 15. Handling performance of semi-active suspensions controlled by different algorithms: (a) time-domain responses; (b) RMS and PP values. 


is almost the same. Since SRF is smaller than URF, the comfort performance is mainly affected by low-frequency responses, 
while the handling performance is mainly affected by high-frequency responses. Therefore, we can see that the LDPWA-Hoo 
controlled suspension performs better RMS and PP results in Fig. 14, while the indicators of the two H^ controllers in Fig. 15 
are basically the same. These phenomena are consistent with the results of the previous conservatism analysis. 

In Fig. 13(a), the SH-ADD controller exhibits the best comfort performance in MFB and HFB, while the LDPWA-Hoo con¬ 
troller performs optimally in the LFB. However, in Fig. 13(b), the SH-ADD controlled suspension displays significantly poor 
performance in the HFB, which is unsatisfactory for the improvement of the handling performance. The RMS and PP results 
in Fig. 15 also illustrate this shortcoming of the SH-ADD algorithm. Since the SH-ADD algorithm is a comfort-oriented con¬ 
troller, the above results are understandable. Conversely, it can be seen from the results in Figs. 13 and 15 that the proposed 
LDPWA-Hoo is capable of achieving good handling performance. 

In summary, the proposed LDPWA-Hoo algorithm shows better comfort performance than the PWA-H^ controller, and at 
the same time, it exhibits better handling performance than the SH-ADD controller. Thus, the multi-objective optimization 
can be achieved. Additionally, once the human sensibility is considered, the comfort performance in the MFB can be further 
improved by increasing the weight parameter p v but this will inevitably cause deterioration of the handling performance in 
the HFB. 

It is worth noting that the RMS indicator and the PP indicator of the SH-ADD controlled suspension in Fig. 14(b) is not 
consistent. Specifically, the RMS value obtained by the SH-ADD controller is obviously better than those of other controllers, 
which can be interpreted according to the PSD results in 13(a). However, the PP indicator obtained by the SH-ADD controller 
is very poor, which can also be found from the time-domain responses in Fig. 14(a). This phenomenon will be further ana¬ 
lyzed in the following simulations. 


4.4.2. Bumpy road simulations 

To analyze the time-domain responses and the control effect under impact disturbance, the simulation condition that the 
vehicle travels through an isolated bump is considered. The road profile is simulated by the following expression: 


Zr(t) 


f(l-cospp)), 

0 , 


ify < f <¥ 

otherwise 


(38) 


where A and L are the height and length of the bump, respectively; l is the length of the flat road ahead of the bump; and V is 
the vehicle speed. These parameters are chosen as that V = 36km/h, A = 0.3m, L = 10 m, and 1= 10 m. Fig. 16 shows the 
time-domain responses and the PP values of signals z s and z u - z r . Correspondingly, the control output of corresponding algo¬ 
rithms is given in Fig. 17. 

As can be seen from Fig. 16, all three algorithms achieve short settling time when the vehicle has passed through the 
bump. However, the PP values obtained by each algorithm are different. Specifically, the PP value resulted from the 
LDPWA-Hoo controller is the smallest, while the peaks obtained by the SH-ADD controller are extremely poor at 1.3s and 
1.7s. In fact, the PP value of the SH-ADD controlled suspension is similar to that of the soft suspension, while the 
PWA-Hoo controlled suspension is closer to the hard suspension. 

The above results can be explained by analyzing Fig. 17. Firstly, inspired by the switching principle of the SH-ADD algo¬ 
rithm, we divide the damping force in the range of Is to 2s into four regions. It can be seen that the damping force output by 
closed-loop algorithms in each region is very close to that of the soft or hard suspension. For example, from Is to 1.25s, all 
the three algorithms behave like the hard suspension, and from 1.25s to 1.5s, their damping forces are closer to that of the 
soft suspension. That is to say, PWA-H^ and LDPWA-H^ controllers also imply a switching mechanism. Further, it can be 
found that each switching behavior in Fig. 17 corresponds to a peak in Fig. 16. To analyze the relationship between the peaks 
and the switching, the damping switching at 1.5s (hard to soft) is taken as an example in Fig. 17. Then, we can see that, for the 
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Fig. 16. Time-domain responses of different control algorithms under an isolated bump. 
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Fig. 17. Damping forces of different control algorithms under an isolated bump: a. damping force in time domain; b. V/F property of PWA-H^ algorithm; c. 
V/F property of SH-ADD algorithm; d. V/F property of LDPWA-H^ algorithm. 


SH-ADD algorithm, this switch happens instantly after a short-time oscillation. However, this oscillation will give rise to a 
conspicuous peak in Fig. 12. For the PWA-H^ algorithm, the switch happens with a terraced reduction in the damping force, 
which is reasonable since the damper performance is restrained by the GQLF. This terraced switch is too slow, which makes 
the force trajectory similar to that of the hard suspension in a long-time range. Unlike the two algorithms, the proposed 
LDPWA-Hoo algorithm carries out this switch with a slope on the damping force, thus the oscillation will not occur, and 
the switch will be faster than that of the PWA-H^ algorithm. Additionally, the damping oscillation here can also explain 
the poor PP performance of the SH-ADD algorithm in Fig. 14. Since the damping oscillation happens with a pure high fre¬ 
quency, the suspension performance reflected by PSD and RMS results in Figs. 13 and 14 will not be influenced. 

In this subsection, the advantage of taking up more available region of the damping force has been verified again and the 
superiorities of the proposed LDPWA-Hoo algorithm have been displayed. However, part of the damping force in Fig. 17 out¬ 
putted by the LDPWA-Hoo algorithm has exceeded the constraint boundary. This is because that the energy of this impact 
disturbance is larger than the design value of w max . A road adaptive controller may achieve better performance, which is 
out of the scope of this paper. 

5. Experiment results 

An experiment system shown in Fig. 18 which includes a test rig (TR) and a measurement and control system (MCS) will 
be used to test the proposed algorithm. The TR can simulate the vertical vibration of suspension systems when a vehicle with 
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Fig. 18. Experiment system for MR semiactive suspensions. 


nominal load is running on sinusoidal road profiles (limited by the ability of this test rig, only frequencies under 6 Hz can be 
tested).Three deflection sensors are assembled in the TR to measure the suspension deflection (y 1 = z s - z u ), the unsprung 
mass displacement (y 2 = z u ), and the road height (simulated by a hydraulic actuator, y 3 = z r ). An accelerometer is used to 
measure the sprung mass acceleration (y 4 = z s ). Parameters of this test rig are same with those used in the previous simu¬ 
lations, and the property of the MR damper coincides with the values given in Table 1 (The passive damping c s is simulated 
by imposing a small shift to the driving current [26]). Equipments which make up the MCS are shown in Fig. 18, where a 
concentrator provides the interaction interface between the TR and the MCS, and a data acquisition (DAQ.) device from 
the IPETRONIK company is used to record the signals which are transmitted by a CAN (Controller Area Network) transceiver 
from the TR and the designed ECU. The power supply of this MCS is a direct current (DC) of 12 V, which can be easily 
obtained from the vehicle battery. Moreover, the ECU satisfies the application requirements of actual vehicles. Specifically, 
a 32-bit MCU (micro-controller unit: MPC5643L, 120 MHz) from the NXP semiconductors, which addresses the ISO 26262 
safety standard, is employed to run the control algorithms and a power switch from the ST microelectronics is used to drive 
the MR damper. The power switch is driven by PWM (pulse width modulation) signals with an adjustable duty cycle which 
corresponds to different electric current values. The signals are obtained from the sensors in analog form. Therefore, AD 
(Analog to Digital) modules in the MCU will be used. Further, by defining y = [z s -z u z u z r z s ] T , the feedback signals of the 
Hoc controllers were obtained from the digital signals by using the following filters: 

1 1 / s s + £ . s s + £ 

X) =- T Vi, x 2 =- T y 2 -y 3 ), x 3 =-V(yi +y 2 , * 4 =-Vy 2 

TS + l 1 TS+1 V2 TS + 1 V1 2J TS + 1 2 

where t = 0.015, £ = 0.01,1 / (ts + 1) is a low pass filter for eliminating the high-frequency noise, and (s + £)/(ts + 1) is a dif¬ 
ferential filter. For the SH-ADD control purpose, the sprung mass acceleration is also filtered as that y 4 /(rs + 1). Considering 
the fast response of MR dampers, the system control period of the ECU is set as 1 ms for all the tested algorithms. 

Based on this experiment system, the responses under sinusoidal excitations of different frequencies can be obtained. 
Further, to give a global view on the concerned frequency range, the experimental results are summarized as the fre¬ 
quency-response (FR) characteristic of the suspension system. According to 44], the FR characteristic of nonlinear systems 
can be approximated with a describing function method. Thus, when a nonlinear system is driving by a sinusoidal excitation 
of a pure frequency, the approximated FR magnitude under this frequency can be calculated with the variance gain from the 
disturbance to the system output. For example, the FR magnitude of the sprung mass acceleration (defined as A acc ) under a 
frequency co can be obtained as that: 


Acc0'tu)| = 201og 10 


Jo fe(0) 2 ^ 


\fo 


' T %(t)) 2 dt) 


(39) 


Similarly, the FR magnitude of the tire deflection (defined as A ur ) can be obtained. Based on this method, the FR characteristic 
of the simulation system in Fig. 18, and the experimental system in Fig. 18 can be both obtained. Specifically, Figs. 19(a) and 
19(b) are the simulated FR results from z r to z s and z u -z r , respectively, and Figs. 20(a) and 20(b) are the corresponding 
experimental results. 

In Fig. 19, the LDPWA-Hoo algorithm shows better comfort and handling performance than the SH-ADD algorithm in LFB. 
Additionally, the performance of the LDPWA-Hoo algorithm in MFB is obviously better than that of the PWA-Hoo algorithm 
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(a) 


(b) 


Fig. 19. Approximated frequency responses of the simulation system: (a) z s ; (b) z u — z r . 



Frequency (Hz) Frequency (Hz) 


(a) (b) 

Fig. 20. Approximated frequency responses of the experiment system: (a) z s ; (b) z u -z r . 

and can approach the effect of the SH-ADD controlled suspension. These results are consistent with the those obtained 
through the random road simulations. More importantly, the same conclusion can also be derived from experimental results 
in Figs. 20. In HFB, Fig. 19(b) shows that the LDPWA-Hoo and PWA-H^ algorithms can both attenuate the vibration of 
unsprung mass, while the performance of SH-ADD algorithm is unsatisfactory for the enhancement of handling performance. 
Therefore, a cross-validation has been given for the proposed algorithm under the whole concerned frequency range. Addi¬ 
tionally, the experiment results obtained through an actual ECU also confirm the applicability of the proposed LDPWA-Hoo 
algorithm for a practical application environment. 


6. Conclusion 

An LDPWA-Hoo algorithm for MR semi-active suspensions has been given in this paper to deal with the load variation and 
road disturbance. To guarantee the designed system performance, the nonlinear constraint on the damping force is handled 
with a piecewise approximation method and a piecewise affine controller. Further, the PDPQLF is used to reduce the system 
conservatism and to obtain a load-adaptive property for the closed-loop system. Simulation results under random and 
bumpy road profiles have been studied. Specifically, the proposed algorithm can reveal a stiff damping property for overload 
conditions and a soft damping property for underload conditions. Further, compared with the GQLF-based PWA-H^ algo¬ 
rithm, the proposed algorithm can utilize the available force region of MR dampers much better, which results to better com¬ 
fort and handling performance. While comparing with the classic SH-ADD algorithm, the proposed algorithm can achieve 
better handling performance under random road profiles and avoid the oscillation of damping force under bumpy roads. 
In the end, the superiority of the proposed algorithm under LFB, which results from the full utilization of MR dampers, 
has been verified by experiments. 
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Appendix A 


Lemma 1. [46]: Let 1(P) be a matrix expression that may relate to P, or may have nothing to do with P. Then, there exists a 
positive-definite matrix P > 0 such that 


--a(p) 

PA 


★ 

P 


<0 


if and only if there exist a positive-definite matrix P > 0 and a matrix R such that 

r-* ( p) * i <0 

r t a p - r - r t < 
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